function [ysgolay] = ex2sgolay()
    load('Data/simdata.mat','t','y','ynf','ysin');
    
    ysgolay = sgolayfilt(y,3,41);   % Apply 3rd-order filter

    figure
    plot(t,ynf,'r',t,ynf+ysin,'m',t,ysgolay,'b');
    legend('ynf', 'ynf+ysin', 'ysgolay')
    title('Sgolay');
    
    save('Data/sgolay.mat','ysgolay');
end